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We propose an approach to realize an n-qubit controlled-[/ gate with superconducting quantum 
interference devices (SQUIDs) in cavity QED. In this approach, the two lowest levels of a SQUID 
represent the two logical states of a qubit while a higher-energy intermediate level serves the gate 
manipulation. Our method operates essentially by creating a single photon through one of the 
control SQUIDs and then performing an arbitrary unitrary transformation on the target SQUID 
with the assistance of the cavity photon. In addition, we show that the method can be applied to 
■ implement an n-qubit controUed-f/ gate with atomic qubits in cavity QED. 

^ I PACS numbers; 03.67.Lx, 85.25.Dq, 42.50.Dv 
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X) ■ I. INTRODUCTION 

Superconducting devices such as Cooper pair boxes, Josephson junctions, and superconducting quantum interference 
devices (SQUIDs) have attracted much attention in the quantum-information community. Because they are relatively 
easy to scale up and have been demonstrated to have relatively long decohcrence times [1-7], they have been considered 
. as promising candidates for physical implementation of quantum computing. It is known that the building blocks of 
quantum computers are single-qubit logic gates and two-qubit logic gates [8]. In the past few years, for SQUID systems, 
, many methods for realizing a single-qubit arbitrary rotation gate and a two-qubit controUed-NOT (or controUed-phase 
C\| shift) gate have been presented [9-16]. 

Multiqubit controlled gates play an important role in constructing quantum computational networks and realizing 
quantum error correction protocols and implementing quantum algorithms. Recently, much attention is paid to physi- 
cal realization of multiqubit controlled gates [17-19]. It is known that when using the conventional gate-decomposition 
protocols to construct a multiqubit controlled gate [20-22], the procedure usually becomes complicated as the number 
of qubits increases. Therefore, it is important to find a more efficient way to implement multiqubit controlled gates. 

In this paper we focus on how to realize an n-qubit controlled-?/ gate with n SQUIDs (l,2,...,n) based on cavity 
QED. Recently, it has been predicted that the strong coupling limit of cavity QED, which is difficult to achieve with 
atoms in a microwave cavity, can readily be realized with superconducting charge qubits [23,24], superconducting 
^ , flux qubits [25], or semiconducting quantum dots [26]. And more recently, the strong coupling cavity QED has been 
C7^' experimentally demonstrated with superconducting charge qubits and flux qubits [27,28] and semiconductor quantum 
dots embedded in a microcavity [29-31]. The controlled gate considered in this paper is shown in Fig. 1, which is 
defined as follows: (a) It contains n—1 control qubits and one target qubit, (b) It leaves the state of the target qubit 
unchanged if not all the control qubits are in the state ]1). (c) However, when all control qubits are in the state ]1) , 
' an arbitrary unitary transformation U is performed on the target qubit. 

To implement the general multiqubit controlled gate described above, three levels ]0) , ]1) , and ]2) of each SQUID 
will be employed (Fig. 2). For each SQUID, the two lowest levels JO) and Jl) represent two logical states of a qubit 
while the higher-energy level ]2) is used to facilitate coherent control and manipulation of quantum states of the qubit. 
The method presented here operates essentially by (a) creating a single photon through one of the control SQUIDs, 
(b) performing a general U on the target SQUID with the aid of the cavity photon as follows [32] 

U = e'^R, (13) Ry (7) R, (S) , (1) 

and (c) finally performing operations to have the cavity mode return to its original vacuum state. In equation (1), 
a,/3, 7, and S are arbitrary real numbers, Ry and R^ represent rotations along y and z axes on a Bloch sphere, which 
are described by matrices 

R.W = [ Q J/2 ) ,^ = P.S, (3) 
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in a single-qubit computational subspace formed by the two logical states |0) = (1,0)"^ and |1) = (0, 1)^ of the target 
qubit. 

As shown below, this scheme has the following advantages: (i) No auxiliary SQUIDs or measurement is needed 

during the entire operation, thus the hardware resources is reduced and the operation is simplified; (ii) As tunneling 
between the qubit levels |0) and 1 1) is not required during the operation, decay from the level 1 1) can be made negligibly 
small during the operation (via prior adjustment of the potential barrier between the qubit levels |0) and |1) [33]) and 
therefore the storage time of each qubit can be made much longer; (iii) The coupling constants of SQUIDs with the 
cavity mode could be different, hence neither identical SQUIDs nor exact placement of SQUIDs is needed; and (v) 
More importantly, the gate operations are significantly simplified as the number of qubits increases, when compared 
with the conventional gate-decomposition protocols. In addition, it is interesting to note that the method can readily 
be extended to obtain an n-qubit controlled-[/ with atomic qubits in cavity QED. 

This paper is outlined as follows. In Sec. II, we review the basic theory of a SQUID coupled to a single-mode cavity 
or driven by a classical microwave pulse. In Sec. Ill, we show a way to realize a two-qubit controUed-t/ gate with two 
SQUIDs coupled to a cavity. In Sec. IV, we discuss how to extend the method to achieve an n-qubit controlled-?/ 
gate with n SQUIDs in cavity QED. In Sec. V, we compare the present method with conventional gate-construction 
protocols. In Sec. VI, we give a brief discussion on experimental issues for the realization of an n-qubit controUed- 
rotation gate. In Sec. VII, we further show how to apply our method to implementing an n-qubit controlled-C/ gate 
with n atoms using one cavity only. A concluding summary is given in the last section. 

II. BASIC THEORY 



The SQUIDs throughout this paper are rf SQUIDs each consisting of a Josephson tunnel junction in a supercon- 
ducting loop (typical size of an rf SQUID is on the order of 10 — 100 /um). The Hamiltonian of an rf SQUID (with 
junction capacitance C and loop inductance L) has the usual form [5] 

Hs = ^+^-^^=;?^-Ejcos(2ni-), (4) 
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where $ is the magnetic flux threading the ring, Q is the total charge on the capacitor, $x is the external magnetic 
flux applied to the ring, and Ej = Ic^o/2tt is the maximum Josephson coupling energy (/c is the critical current of 
the junction and = h/2e is the flux quantum). 

A. SQUID-cavity resonant interaction 

Consider a SQUID coupled to a singlc-modc microwave cavity field. The SQUID is biased properly to have a 
A-type configuration formed by three lowest levels, denoted by |0) , |1) and |2) with energy eigenvalues Eq, Ei, and 
E2, respectively (Fig. 2). The transition frequency between the two levels \i) and \j) is Uij = cuij/ (27r) = \Ei — Ej\ /h 
€ {0, 1,2} ,z 7^ j). Suppose that the coupling of |0) , |1) and |2) with other levels of the SQUID via the cavity 
is negligible, which can readily be achieved by adjusting the level spacings of the SQUID [33]. We can show that 
when the cavity mode is resonant with the |0) |2) transition while decoupled (highly detuned) from the |1) |2) 
transition and the |0) ^ |1) transition of the SQUID, the interaction Hamiltonian in the interaction picture, after the 
rotating-wave approximation, is described by [12] 

Hi = n{ga+\0){2\+h.c.), (5) 

where a"*" and a arc the creation and annihilation operators of the cavity mode, and g is the coupling constant between 
the cavity mode and the |0) ^ |2) transition of the SQUID. For a superconducting one-dimensional transmission line 
standing- wave cavity, g is given by 

n9{x) = ^.[^{0\^\2)sm(^-^x), (6) 
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where Mgc is the SQUID-cavity mutual inductance, Lq is the inductance per unit length of the cavity, I is the length 
of the cavity, I'c = (27r) is the frequency of the cavity mode with wavelength A, and x is position of the center of 
the SQUID in the cavity. 

In the case when the cavity is initially in the photon-number state |n), the time evolution of the states of the 
system, under the Hamiltonian (5), is as follows 
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|0) \n) cos ^/ngt |0) \n) — i sin ^/ngt |2) |n — 1) 



\2)\n) -ismy/n+lgt\0) \n + 1) + cosVn + lgt\2) \n) . (7) 

SQUIDs may have non-uniform device parameters and/or be not exactly placed in the cavity. Therefore, the coupling 
strength g may not be identical for different SQUIDs. In the following, we will replace g by gi, g2, and g„ for 
SQUIDs 1,2,..., and n, respectively. 

As shown below, resonant interaction between the cavity mode and the |1) ^ |2) transition of the SQUID is needed. 
In this case, the time evolution for the states of the system is similar to Eq. (7). One just needs to replace |0) in 
Eq. (7) by |1) . The coupling constant of the cavity mode with the |1) ^ |2) transition of the SQUID can be set to 
be the same as that with the |0) ^ |2) transition of the same SQUID by simply changing the external flux bias from 
(0.5 - 5) $0 to (0.5 + 5) $0 (see Fig. 3). 

B. SQUID-cavity off resonant interaction 

Consider a system composed of a SQUID and a single-mode cavity. Suppose that the cavity mode is off resonant 
with the |0) |2) transition (i.e., A = u)2o — u)c » g) while decoupled from the |1) |2) transition and the |0) |1) 
transition of the SQUID [Fig. 4(a)]. Here, A is the detuning between the cavity-mode frequency and the |0) ^ |2) 
transition frequency of the SQUID, and g is the coupling constant of the cavity mode with the |0) ^ |2) transition of 
the SQUID. Under this condition, the SQUID has a negligible probability of making a transition between the ground 
level |0) and the excited level |2) . Therefore, the effective interaction Hamiltonian in the interaction picture can be 
written as [34] 

ife = ?l|(|2)(2|-|0)(0|)a+a. (8) 

From the Hamiltonian (8), it is straightforward to see that if the cavity mode is initially in the photon- number 
state |n), the time evolution of the states of the system is then given by 

|0) \n) ^ e'^"* |0) jn) , 

|2) \n) ^ e-'^"* |2) |n) . (9) 

In the following, we will need to use off-resonant interaction between the cavity mode and the |1) |2) transition 
of the SQUID [Fig. 4(b)]. In this case, the time evolution for the states of the system is similar to Eq. (9). One just 
needs to replace |0) in Eq. (9) by |1) . The detuning is given by A = uj2i ^ ^^c- 9 is the coupling constant between the 
cavity mode and the |1) ^ |2) transition of the SQUID. As described above, i.e., by changing the external flux bias 
from (0.5 — 6) $o to (0.5 + 5) $0) both A and g can be set to be the same as those for the case of the cavity mode 
being off-resonant with the |0) |2) transition of the same SQUID (Fig. 4). 

C. SQUID-microwave resonant interaction 

In this section, we consider a SQUID driven by a classical microwave pulse with the magnetic component B^^(r, t) = 
B^^(r) cos {oj^wt + (!>)■ Here, B^i„(r), 

ijJuwi and are the magnetic field amplitude, carrier frequency, and phase of 
the microwave pulse. Assume that the microwave pulse is resonant with the |1) |2) transition but decoupled from 
the jo) ^ |2) transition and the ]0) |1) transition of the SQUID. The interaction Hamiltonian in the interaction 
picture is then given by 

if, = ^(Qe^^]l)(2l+h.c.), (10) 
where fi is the Rabi frequency of the pulse, which has the following form [12] 

^{t) = ^{l\^\2) jB^^{v)-dS. (11) 
From the Hamiltonian (10), it is easy to find the following state rotation 

|1) ^ cos^t|l) -ze-''^sin^t|2), 

|2) ^ -ie*"^ sin ^t\l) + cos |2) . (12) 
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III. TWO-QUBIT CONTROLLED-C/ GATES 



For two qubits, there are four computational basis states denoted by |00) , |01) , |10) , and |11), respectively. A 
two-qubit controlled- gate is described by 

|00) |00) , 
|01) ^|01), 
|10) -|1)C/|0), 

111) -|1)C/|1), (13) 

which implies that if and only if the control qubit (the first qubit) is in the state |1), a unitary transformation U is 
performed on the target qubit (the second qubit) and nothing happens otherwise. 

Now let us discuss how to obtain the two-qubit controUed-t/ gate (13) with two SQUIDs 1 and 2 coupled to a 
microwave cavity. The SQUIDs considered here have the A-type three-level configuration as depicted in Fig. 2. 
The transition between any two levels for each SQUID is initially decoupled from the cavity mode (e.g., via prior 
adjustment of the level spacings). And the cavity mode is initially in the vacuum state \0)^ . To realize the gate (13), 
it is required to perform the following transformation 

|0)i|0)2®|0), ^|0)J0)2®|0),, 
|0)i|l)2®|0), ^|0)J1)2®|0),, 

10)2 ® |0), ^ |l)i [e'^R. (/3) Ry (7) Rz (S) 10)2] ® |0), , 

|i>i |i>2 ® I0)c ^ |i>i W^'R. W Ry (7) Rz (S) 11)2] ® |0), , (14) 

where subscripts 1 and 2 represent SQUID 1 (the control qubit) and SQUID 2 (the target qubit), respectively. 

We note that the unitary transformations, involved in the last two lines of Eq. (14), can be realized through the 
following operation sequence: 

First, creating a single photon through SQUID 1 as follows 

|l)l|0)e-|0)i|l)e. (15) 

Second, performing rotations Rz{S), Ryi'j), and then Rz{l3) on the states of SQUID 2 with the assistance of the 
photon, i.e., 

o (s). 10)2 |l),-e-^ 10)2 |1),, . . 

11)2 |l).-e^^ 11)2 |1)., 

o . 10)2 |l)c - (cos i 10)2 + sin 5 11)2) |1), , . . 

11)2 |l)e-(-sini 10)2 +cosi 11)2)11),, ^''^ 

o (o. 10)2 |1),- 6-/3 10)2 11),, . . 

^^(^^ |l),|l),^e»^|l),|l),. (1^^ 

Third, performing a phase-shift e'" on the states of SQUID 2 with the aid of the photon, i.e., 

|0)2|1), -e^"|0)2|l),, (19) 
11)2 |1), -e*« 11)2 |l)e. 

Last, returning SQUID 1 and the cavity mode to their original states, i.e., 

|0)i|l)e-|l)i|0),. (20) 

In the following, we will list operations required for the realization of the above transformations (15)-(20). 

Step (i): Apply a n microwave pulse (rir^™ = tt, where r^^ is the pulse duration) with 6 = —n/2 to SQUID 1 [Fig. 
5(a)]. The pulse is resonant with the |1) |2) transition of SQUID 1. After the pulse, the transformation |1) |2) 
of SQUID 1 is obtained. 
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Step (ii): Bring the |0) ^ |2) transition of SQUID 1 to resonance with the cavity mode for an interaction time 
n = T^l (251) [Fig. 5(b)], resulting in 12), |0)^ ^ -i |0), 

After the operations of Step (i) and Step (ii), the transformation (15) is obtained as follows 

|l>i|0),^|2)jO),^^-»|0)Jl), (21) 

up to a phase factor — i, which is inevitable according to Eq. (7) but can be removed by introducing a phase factor i 
to the transformation (20) (see below). 

Step (iii): Adjust the level structure of SQUID 2 and apply a tt microwave pulse with <^ = — 7r/2 to SQUID 2 [Fig. 
5(c)]. The pulse is resonant with the |1) |2) transition of SQUID 1, leading to the transformation |1) — > |2) of 
SQUID 2. 

Step (iv): Adjust the level structure of SQUID 2 to obtain off-resonant interaction between the cavity mode and 
the jo) ^ |2) transition of SQUID 2 [Fig. 5(d)]. After an interaction time ti = SA/g^, the state [0)2 [1)^ goes to 
e-'^ [0)2 |1)^ while the state [2)2 |1)^ changes to e'^ \2)^ |1)^ . 

Step (v): Repeat the operation of Step (iii) but set (f) = 7r/2, leading to the transformation |2) — > |1) of SQUID 2. 

It is easy to see that after the operations of Steps (iii)-(v), the transformation (16) is implemented as follows: 

10) 2 10)2 |1), '-^^ e-^' 10)2 |1), e-^' 10)2 |1), , 

11) 2 12)2 |1), e^' 12)2 |1>, e^' 11)2 |1>, , (22) 

Step (vi) Bring the |1) ^ |2) transition of SQUID 2 to resonance with the cavity mode for an interaction time 
T3 = tt/ (2(72) [Fig- 5(e)]. As a result, the state [0)2 remains unchanged while the state [1)2 [1)^ changes to 
-i|2)2|0),. 

Step (vii): Bring the |0) ^ |2) transition of SQUID 2 to resonance with the cavity mode for an interaction time 
Ti = 7/(2,92) [Fig. 5(f)], resulting in [0)2 ^ cos ||0)2 |1), - i sin ^ [2)2 |0)^ and [2)2 |0)^ ^ sin ||0)2 |1), + 
cos ^12)210),^. 

Step (viii): Bring the |1) |2) transition of SQUID 2 to resonance with the cavity mode for an interaction time 
T5 = Stt/ (25(2) [Fig. 5(e)]. As a result, the states [0)2 [1)^ remains unchanged while the state [2)2 |0)^ becomes 

*|l)2|l)c- 

One can see that after the operations of Steps (vi)-(viii), the transformation (17) is obtained as follows: 

10) 2 |l)c - 10)2 |l)c cos I 10)2 |1), - i Sin 1 12)2 |0), ^ (cos I 10)2 + sin | [1)2) |1), , 

11) 2 |1), -z 12)2 |0), - sin 1 10)2 |1), - i cos 1 12)2 |0), (- sin 1 \0), + cos | [1)2) |1), • 

(23) 

Steps (ix)-(xi): Repeat the operations of Steps (iii)-(v) but set the cavity-SQUID off-resonant interaction time as 

tq = f3A/g^, leading to the transformation (18). 

Step (xii): Adjust the level structure of SQUID 2 to obtain off-resonant interaction between the cavity mode and the 
|0) |2) transition of SQUID 2 [Fig. 5(g)]. After an interaction time T7 = the state [0)2 11)^ goes to e*" ]0)2 

while the state [1)2 [1)^ remains unchanged. 

Step (xiii): Adjust the level structure of SQUID 2 to achieve off-resonant interaction between the cavity mode and 
the |1) ^ |2) transition of SQUID 2 [Fig. 5(h)]. After an interaction time rg = the state [1)2 [1)^ changes to 

e*" [1)2 ll)(, while nothing happens to the state [0)2 ll)^ • 

It can be seen that after the operations of Step (xii) and Step (xiii), the transformation (19) is obtained as follows: 

10)211),*^^ e»l0)2ll),''i'^ 10)2 11),, 

|l>2|l)c*^V>2|l>c''-^'^e-ll)2ll),. (24) 

Step (x-iv) Bring the ]0) ^ \2) transition of SQUID 1 to resonance with the cavity mode for an interaction time 
rg = Stt/ (251) [Fig. 5(b)], resulting in \0)^ ^ i \2)^ [0), . 
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Step (x-v): Apply a tt microwave pulse with cf) = tt/2 to SQUID 1 [Fig. 5(a)]. The pulse is resonant with the 
|1) |2) transition of SQUID 1. After the pulse, the transformation |2) |1) of SQUID 1 is obtained. 
One can see that the operations of Step (x-iv) and Step (x-v) lead to the following transformation 

|0)i|l),^z|2)j0),^z|l),|0),, (25) 

which is actually the transformation (20) up to a phase factor i. 

In above, we have explicitly shown how to realize the transformations (15)-(20), i.e., performing a general U on 
SQUID 2 (the target) when SQUID 1 (the control) is initially in the state |1) . On the other hand, it is noted that 
the following states of the system 

|0)i 10)2 |0),,|0), 11)2 |0), (26) 

remain unchanged during the entire operation. This is because: (a) During the operation of Step (i), the state |0) of 
SQUID 1 was not affected by the applied microwave pulse, since the |0) |2) transition and the |0) ^ |1) transition 
of SQUID 1 are decoupled from the pulse; and (b) No photon was emitted to the cavity during the operation of Step 
(ii), when SQUID 1 is initially in the state |0). Hence, it can be concluded that the transformation (14), i.e, the 
two-qubit controlled-[/ gate (13) was implemented with two SQUIDs after the above manipulation. 

Before closing this section, several issues need to be addressed. The irrelevant SQUIDs in each step of the operation 
need to be decoupled from the cavity/pulse during the cavity/pulse-SQUID interaction. The cavity mode needs to 
be not excited during the SQUID-pulse resonant interaction. In addition, for each SQUID, the coupling of the levels 
|0), |1) , and |2) with the other levels should be negligible. Note that for a SQUID, the level spacings can readily be 
changed by varying the external flux or the critical current Ic (e.g., for variable barrier rf SQUIDs) [33]. Therefore, 
these conditions can in principle be satisfied by adjusting the level spacings of the SQUIDs. 

Imperfect decoupling between the irrelevant SQUIDs and the cavity during the operations using off-resonant in- 
teraction could, in principle, result in gate errors. Note that the population of the level \j) of any irrelevant SQUID 
initially in the state induced due to the coupling between the \i) \j) transition and the cavity mode, is on the 
order of Pj ~ gfj/ {gfj + Af^ ), where i,j G {0, 1, 2} and i ^ j, gij is the coupling constant of the cavity mode with 
the \i) ^ \j) transition of the SQUID, and Ajj = LUij — ujc is the detuning of the \i) ^ \j) transition from the cavity 
mode. Therefore, we remark that the coupling between the i£relevant SQUIDs and the cavity can be made negligible 
as long as the condition Aij >> gij and Aij/gfj » A/g^, A/g^ can be satisfied. Here, "g and g are the off-resonant 
coupling constants described above, i.e.. the coupling constants of the cavity mode with the transition between the 
corresponding two levels of the SQUID which is involved during the operation using off-resonant interaction. The 
required decoupling condition can be obtained by adjusting the level spacings of the irrelevant SQUIDs before the 
off-resonant operations, so that the transition frequency between any two levels of the irrelevant SQUIDs is highly 
detuned from the cavity-mode frequency. It can be achieved with the available experiment technique because the 
level spacings of a SQUID can be adjusted rapidly in experiment (~lns). Since a more quantitative answer to the 
question of "how well decoupled the irrelevant SQUIDs need to be from the cavity" requires a very lengthy and 
complex analysis, we will not give a detailed discussion. 

IV. N-QUBIT CONTROLLED-;/ GATES 

For n qubits, there are a total number of 2" computational basis states from [00. ..0) to [11. ..1), which form a set 
of complete orthogonal bases in a 2"-dimensional Hilbert space of the n qubits. As discussed in the introduction, an 
n-qubit controlled- /7 gate performs an arbitrary unitary transformation U on the target qubit only when the n — 1 
control qubits are all in the state ]1), i.e., 

|o) ^ |l)®("-i)C7]0), 

^ (27) 

while nothing happens to all other 2" — 2 computational basis states. In Eq. (27), the first n — 1 qubits represent 
control qubits while the last qubit acts as a target. In the following, we will discuss how this gate can be achieved 
with n SQUIDs coupled to a cavity. 

The n SQUIDs are labeled by 1, 2, and n. The first n — 1 SQUIDs (1, 2, n — 1) represent control qubits while 
SQUID n is the target qubit. Suppose that all SQUIDs (1, 2, n) are initially decoupled from the cavity (which is in 
the vacuum state). We find that the n-qubit controlled-C/ gate described above can be obtained through the following 
sequence of operations (from right to left) 
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(28) 



where HILi = f7(„_i)c • • • U2cUic\ Ui denotes the operation on SQUID 1 represented by matrix 




in the basis states = (0,1)^ and |2)-^ = (1,0)^; Uic is a joint operator on the SQUID / and the cavity mode 
(Z = 1,2, ...,n — 1), represented by the matrix 





which performs an arbitrary unitary transformation U on the SQUID n while nothing to the cavity state. 
From the description in the previous section, it can be seen that; 

(i) U\ (Ui) can be reahzed by a tt microwave pulse {flr^w = tt, where Tij,yj is the pulse duration) with (f) = —■jt/2 
{tt/2) and o;^^ = 0^21 to SQUID 1; 

(ii) Uic corresponds to the operation of bringing the |0) ^ |2) transition of SQUID I {I = l,2...,n— 1) to resonance 
with the cavity mode for an interaction time Ti = w/ {2gi); 

(iii) Unc can be realized via the operations of Steps (iii)-(xiii) described in the previous section. For the present 
case, the qubit involved in the operations is SQUID n instead of SQUID 2. 

(ii) corresponds to the operation of bringing the |0) |2) transition of SQUID I {I = 1,2..., n— 1) to resonance 
with the cavity mode for an interaction time rj = Stt/ (2gi) . 

To further understand Eq. (28), let us give some explanation on system evolutions after the operations described 
above. 

First, the state |1) of SQUID 1 is changed to the ground state |0) and a single photon is emitted to the cavity mode 
after the joint operation UicUi. 

Second, after the joint operation nr=2^ ^'c, the photon is left in the cavity mode when SQUIDs (2, 3, n — 1) are 
initially in the computational basis state |ll...l) ; however, it is absorbed by SQUIDs (2,3, ...,n — 1) initially in all 
other computational basis states. 

Third, after the operation Unc, an arbitrary unitary transformation U is performed on SQUID n with the assistance 
of the cavity photon; in contrast, nothing happens to the SQUID n if no photon, after the previous operation Y\^Z2 ^ic, 
is left in the cavity mode. 

Fourth, after the joint operation 111^2^ ' photon, originally absorbed by SQUIDs (2, 3, n — 1), is emitted 

back to the cavity mode. 

Last, after the joint operation U^U^^, the cavity mode returns to its original vacuum state \0)^. and the state |0) of 
SQUID 1 is changed back to the initial state |1) . 

It should be mentioned that when SQUID 1 is initially in the state |0), nothing happens to the whole system during 
the entire operation, due to the reason discussed in the previous section. 

To see it more clearly, let us consider the case of three qubits, i.e., realizing a three-qubit controlled-?/ gate with 
three SQUIDs (1,2,3) in cavity QED. For three qubits, there are a total number of eight (2'') computational basis 
states, denoted by |000) , |001) , jlll) , respectively The first four basis states |000) , |001) , jOlO) , and |011) of the 
three SQUIDs remain unchanged during the operations described above. However, if the three SQUIDs are initially 
in the other four basis states, the states of the whole system after each unitary transformation of Eq. (28) (n = 3) 
are: 



|100)|0), |200)|0), 
|101)|0), |201)|0), 



i|000)|l), 
«|001)|1), 



- 1020) |0), 
- 1021) |0), 
*|010)|1), 
*|011)|l)e 



|110)|0), |210)|0), 
|111)|0), |211)|0), 



i|010)|l), 
i|011)|l)e 
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-|020)|0), -*|000)|1), |200)|0), |100)|0), 

-|021)|0), ut -x|001)|l), ut |201)|0), ut |101> |0), 

101) {U |0)) |1), ^ 101) iU |0)) |1), |21) {U |0)) |0), |11) {U |0)) |0), ' ^-"^^ 

-i |01) (t/ |1)) |1), -i |01) (C/ |1)) |1), |21) ([/ |1)) |0), 111) ([/ |1)) |0), 

where |«jfc) is abbreviation of the state \j)2 |fc)g of SQUIDs (1,2,3) with i,j,k G {0, 1,2}. From Eq. (32), it can 
be seen that a three-qubit controlled-t/ gate is achieved with three SQUIDs, where the third SQUID is the target 
qubit, after the transformations described in Eq. (28). 

On a final note, we point out that a single-mode cavity is not necessary since for a multi-mode cavity one can in 
principle choose one mode to interact with the SQUIDs while have all other modes well decoupled from the three 
lowest levels of the SQUIDs. In addition, the method presented here is applicable to a ID, 2D, or 3D microwave 
resonator /cavity as long as the conditions described above are satisfied. 

V. COMPARISON WITH CONVENTIONAL GATE CONSTRUCTION 

The universality of quantum computation implies that it is possible to generate arbitrary n-qubit gates by using 

sequences of onc-qubit and two-qubit gates only [35-37]. Barenco et. al. have first developed methods to design 
networks for n-qubit controlled gates [20]. They have shown that it requires 2""^ — 1 two-qubit controlled-^ and 
controlled- F"*" gates and 2"~^— 2 two-qubit controlled-NOT gates to accomplish an n-qubit controlled-?/ gate described 
above (n > 3). Here, V satisfies V'^ = U. Namely, at least 2" — 3 steps of operations are required, assuming that 
realizing any two-qubit controlled gate requires one-step operation only. 

Since the work of Barenco et. al., much attention has been paid to optimal implementation of quantum gates. 
Recently, Mottonen et. al. have considered a generic elementary gate sequence for realizing a general multiqubit gate 
[21]. More recently, Bergholm et. al. have presented quantum circuits with uniformly controlled one-qubit gates [22]. 
According to their results, 2"~^ — 1 two-qubit controlled-NOT gates, 2"~^ one-qubit gates, and a single diagonal 
n-qubit gate arc needed to construct the above n-qubit controUed-C/ gate. Hence, at least 2" steps of operations are 
required, provided that any two-qubit gate, one-qubit gate, or an n-qubit diagonal gate can be implemented using 
one-step operation only. 

It is interesting to note that when compared with the use of the conventional gate-decomposition protocols [20,22], 
our method significantly reduces the number of operations needed to implement the above n-qubit controlled-?/ gate. 
As shown above, our method only needs 2n-|-ll steps of operations. The advantage of this method starts at n = 5 and 
becomes more dramatic as n increases. It should be mentioned that the comparison presented here is a conservative 
case because it was assumed above that a two-qubit controlled-F/ y+/NOT or a n-qubit diagonal gate can be realized 
with one-step operation only. 

VI. DISCUSSION 

In this section, we give a brief discussion on relevant experimental issues for the implementation of n-qubit controUed- 
U gates. Without loss of generality, let us consider n identical SQUIDs (1,2, ...,n) at locations where the Be fields 
arc the same (e.g., antinodes of the cavity field). Thus, we have gi = g {I = 1,2, ...,n). The total operation time is 

given by 

r = [2n + 7/ (27r)] t^^^ + 2a4^^ + {(3 + 5)4^^ + (2n + 9)r„ + 4r^^, (33) 

where ri^^ = 'K/g,T^'' = A/g"^, ri^^ = A/g^, and is the typical time required for adjusting the level spacings of a 
single SQUID. For the method to work, r should be much shorter than the energy relaxation time 7^^^ of the level |2) , 
and the lifetime of the cavity mode k,~^ = Q/2tti'c, where Q is the (loaded) quality factor of the cavity. In addition, 
direct coupling between SQUIDs needs to be negligible since this interaction is not intended. 

These requirements can in principle be realized, since one can: (i) reduce ri^'' by increasing the coupling constant 
g, (ii) shorten by rapid adjustment of the level spacings of the SQUIDs, (iii) increase by employing a high-Q 
cavity so that the cavity dissipation is negligible during the operation, and (iv) design SQUIDs and control (readout) 
circuitry so that the energy relaxation time the level |2) is sufiiciently long. In addition, it is noted that direct 

interaction between SQUIDs can be made negligible as long as the following condition is satisfied 

M,, < Msc, (34) 

where Mgs is the mutual inductance between two adjacent SQUIDs and Mgc is the mutual inductance between each 
SQUID and the cavity. 
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For the sake of definitivcness, let us consider the experimental feasibihty of reahzing a five-qubit controUed-/7 
gate using SQUIDs with the parameters hsted in Table 1. A five-qubit controUed-t/ gate performs the following 
transformation 

11111) |0) |llll)C/(a,/?,7,^)|0), 

^ 11111) C/(a, /3, 7, <5) 1 1), (35) 

when the four control qubits are in the state |1111) while nothing otherwise. Here, 7/2 e [0,27r], and a,(i/2,5/2 G 
[— TT, tt] (The positive or negative value for a,/?, and S can be achieved by setting blue or red detuning). Note that 
SQUIDs with the parameters in Table 1 are readily available at the present time [2,3,38] and have the desired three-level 
structure as depicted in Fig. 2. For a superconducting one dimensional standing-wave CPW (coplanar waveguide) 
cavity with the parameters listed in Table 1 and SQUIDs placed along the cavity axis (Fig. 6), one has Mgc ^ 10^ 
pH. When each SQUID is located at one of the antinodes of the cavity mode (Fig. 6), a simple calculation gives 
g 5.8x 10^ s~^, resulting in ri^^ 0.5 ns. On the other hand, as a rough estimate, we assume 5 ~ g ^ O.Sg, A ~ lOg, 
and A ~ lOg, which can be readily achieved by adjusting the level spacings. As a result, we have Tc ^ Tc ~ 3.4 
ns. With the choice of r^tu ^ Ta ^ Tc^\ one has r ^ 81.1 ns for 7/2 = 27r and a, 13/2, 7/2, S/2 ~ ±7r (a case requiring 
the longest operation time), which is much shorter than 72"^ ~ 3.2 /xs and ~ 0.8 fis for a cavity with Q ~ 6 x 10"^. 
Note that superconducting CPW resonators with a quality factor of Q > 10^, patterned into a thin superconducting 
film deposited on the surface of a silicon chip, has been experimentally demonstrated [39] (also see Refs. [23,26,27] 
regarding its application for loaded superconducting qubits or semiconductor qubits). 

For a cavity with Vc = 11.4 GHz, the wavelength of the cavity mode is A ^ 10.5 mm. For each SQUID being placed 
at an antinode of the Be field (Fig. 6), one has D ~ 5.25 mm, where D is the distance between the two nearest 
SQUIDs. A simple numerical calculation gives Mgs ~ 0.1 aH, which is much smaller than M^c- Hence, the condition 
of negligible direct coupling between SQUIDs is very well satisfied. 

Our above analysis demonstrates that the realization of a five-qubit controUed-C/ gate is possible using SQUIDs 
and a cavity within the present technology. In addition, we point out that a quantum controlled Hadamard gate 
with a larger number of control qubits can in principle be obtained by increasing the length of the cavity though the 
conditions of r <C 7^^, becomes increasingly difficult to satisfy. 

VII. iV-QUBIT CONTROLLED-?/ GATES WITH ATOMS 

In this section, we discuss how to extend the above method to realize an n-qubit controlled-?/ gate with three-level 
atoms, by the use of one cavity. 

Consider n identical atoms (1, 2, ...,n) each having A- type level configuration formed by two ground states and an 
excited state (Fig. 7). In accordance with the previous section, we use |0) and |1) to represent the two ground states 
and |2) to indicate the excited state. The dipole transition between |0) and |1) is forbidden due to the definite parity 
of the wave function. In the following, the two logical states of a qubit are represented by the two ground states |0) 
and |1) of each atom. And, atoms (1, 2, n — 1) act as control qubits while atom n is the target qubit. 

We note that an n-qubit quantum controUed-t/ gate described by Eq. (27) can be realized with three-level atoms 
using the prescription of Eq. (28), by performing the following operations: 

(a) . Apply a tt/2 classical pulse with (f> = — 7r/2 to atom 1, resonant with the |1) ^ |2) transition of atom 1. This 
pulse leads to the transformation ]1) — > |2) of atom 1, which accomplishes the transformation Ui of Eq. (28). 

(b) . Send atom 1 through the cavity. The cavity mode is initially in the vacuum state |0)^ and resonant with the 
jo) 1 2) transition of atom 1 [Fig. 7(a)]. Choose the atomic velocity appropriately so that the passage time of atom 
1 through the cavity equals to tt/ {2gi). Thus, after atom 1 exits the cavity, the state [0)^ ]0)^ remain unchanged, 
while the state |2)-^ |0)^ changes to —i \0)-^^ |1)^ . This operation implements the transformation Uic in Eq. (28). 

(c) . Send atoms (2, 3, n — 1) through the cavity one after another in a way that no more than one atoms stay in 
the cavity simultaneously. The cavity mode is resonant with the ]0) \2) transition of each atom [Fig. 7(a)]. Choose 
the atomic velocity appropriately so that the duration of atom I in the cavity equals to tt/ {2gi) {I = 2, 3, ...n — 1). As 
a result, there is no change for the state |0); |0)^ , |1); |0)^ , and |1); jl)^ ; while the state )0); ]1)^ becomes —i ]2); [0)^. 
This process completes the transformation 111^2 of Eq. (28). 

(d) . The transformation Unc in Eq. (30) is performed on atom n and the cavity mode, which is realized as follows: 
(d.l) First, apply a 7r/2 classical pulse with = — 7r/2 to atom n, resonant with the |1) ^ |2) transition of atom 

n and resulting in the transformation ]1) |2) . Second, adjust the cavity frequency [40] to obtain an off-resonant 
interaction between the cavity mode and the |0) ^ |2) transition of atom n [Fig. 7(b)] and then send atom n through 

the cavity. After an interaction time the state ]0)^ |1)^ goes to e~"^ lO),^ |1)^ while the state |2)^ |1)^ changes to 
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e"' |2)^ Last, apply a 7r/2 pulse (with (j) = 7r/2) to atom n, resonant with the |1) ^ |2) transition of atom 1 and 
resulting in the transformation |2) |1) . After these operations, a rotation {5) on the two states |0) and |1) of 
atom n is obtained while the cavity mode remains in one-photon state. 

(d.2) First, adjust the cavity frequency so that the cavity mode is resonant with the |1) ^ |2) transition of atom 
n [Fig. 7(c)]. Then send atom n through the cavity for an interaction time it/ {2g'^) so that the state |0)^|1)^ 
does not change while the state goes to —i |2)^ |0)^ . Second, adjust the cavity frequency so that the cavity 

mode is resonant with |0) <-> |2) transition of atom n [Fig. 7(a)]. Then send atom n back through the cavity 
for an interaction time 7/ (2(7„), leading to the rotation |0)„ |1)^ cos ^ |0)^ |1)^ — i sin^ |2)^ |0)^ and |2)^ |0)^ 
— isin ^ |0)^ |1)^ + cos ^ |2)^ |0)^ . Last, have the cavity mode resonant with the |1) <-> |2) transition of atom n via 
the adjustment of the cavity frequency [Fig. 7(c)]; and then send atom n through the cavity for an interaction time 
Stt/ {2gn) ■ As a result, the state 10),^ |1)^ remains unchanged while the state |2)^ |0)^ becomes i |1)^ |1)^ . It can be 
seen that after these operations, a rotation gate Ry (7) on the two states |0) and |1) of atom n is realized in the same 
manner as shown in Eq. (23). 

(d.3) To achieve a rotation R^ {13) on the two states |0) and |1) of atom n, one just needs to perform the same 

operations as described in (d.l) by simply setting the time of atom n crossing the cavity as 

(d.4) A common phase e*" for the two states |0) and 1 1) of atom n can be obtained as follows. First, adjust the cavity 

frequency to obtain an off-resonant interaction between the cavity mode and the |0) ^ |2) transition of atom n [Fig. 
7(d)] and then send atom n through the cavity for an interaction time S^, resulting in |0)„ |0)„ |1)^ while 

no change for the state ll)^^^ |1)^ . Second, adjust the cavity frequency to obtain an off-resonant interaction between 
the cavity mode and the |1) |2) transition of atom n [Fig. 7(e)]. Then, send atom n back through the cavity for 
an interaction time which leads to ]1)„ e*" ]1)^ . 

The operations described above have accomplished a general transformation U on the target qubit (atom n) with 
the assistance of the cavity photon. However, it is noted that the original states of the cavity mode and atoms 
(1, 2, ...n — I) have also changed after the above operations. Therefore, it is necessary to return the cavity mode and 
atoms (1, 2, n — 1) (control qubits) to their original states, which can be done through the following operations: 

(e) Send atoms (2, 3, n — 1) through the cavity one after another, without more than one atoms staying in the 
cavity simultaneously. The cavity frequency is adjusted to have the cavity mode resonant with the ]0) ^ ]2) transition 
of each atom [Fig. 7(a)] . Choose the atomic velocity appropriately so that the duration of atom I in the cavity equals 
to 3tt/ {2gi) {I = 2, 3, ...n — 1). As a result, nothing happens to the states JO); JO)^ and ]1); ]0)^ while the state 12)^ ]0)^ 
becomes ^10)^ jl)^ {I = 2, 3,...n— 1). After this operation, atoms (2,3, ...n — 1) return to their original states and 
the photon, originally absorbed by atoms (2,3, ...,n — 1), is emitted back to the cavity. This process realizes the 
transformation 1 ^ic '^^ (28)- 

(f) . Send atom 1 back through the cavity. The cavity mode is resonant with the ]0) ]2) transition of atom 1 
[Fig. 7(a)]. Choose the atomic velocity appropriately so that the passage time of atom 1 through the cavity equals 
to 37r/ (2gi) . Thus, after atom 1 exits the cavity, the state jO)-^ )0)^ remains unchanged but the state jO)-^ changes 
to i \2)^ )0)^ . After this operation, the cavity mode returns to the original vacuum state |0)^ and the transformation 
Utc of Eq. (28) is obtained. 

(g) . Finally, apply a 7r/2 classical pulse (with (j) = 7r/2) to atom 1. The pulse is resonant with the )1) )2) 
transition of atom 1. This pulse leads to ]2) |1) , completing the transformation Ui of Eq. (28). 

The present scheme has the following advantages: 

(i) No adjustment of the level spacings for each atom is required during the operations; 

(ii) Only one cavity is required; 

(iii) No identical atom-cavity coupling constants are needed; and 

(iv) The total number of basic operations is 2n + 11, which is much less than that required by the conventional 
gate-decomposition protocols [20,22] when n is a larger number (n > 5). 

VIII. CONCLUSION 

We have presented a method to realize a multiqubit quantum controlled-C/ gate with SQUIDs coupled to a microwave 
cavity. The method operates essentially by creating a single photon through one of the control SQUIDs and then 
exchanging the photon between the control SQUIDs and the cavity mode before and after a unitary transformation U 
is performed on the target SQUID. The method has these advantages: (i) Since no tunneling between the qubit levels 
)0) and Jl) is required, decay from the level Jl) can be made negligibly small during the operation, via prior adjustment 
of the barrier of the double- well potential [33]; (ii) As neither measurement on SQUIDs/photons nor auxiliary SQUID 
is needed, the operation is simplified and hardware resources is saved; (iv) Because coupling constants of SQUIDs with 
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the cavity are not required to be identical, inevitable nonuniformity in device parameters is tolerable and non-exact 
placement of SQUIDs is allowed; (iv) The method can in principle be applied to obtain an n-qubit controlled-[/ gate 
with a large number n, and (v) More interestingly, the gate operations are significantly simplified as the number of 
qubits increases, when compared with the use of the conventional gate-decomposition protocols. As shown above, the 
present method can be extended to implement a multiqubit controUed-C/ gate with atoms in cavity QED. Finally, it 
is noted that the method is also applicable to the realization of a multiqubit controlled-C/ gate with quantum dots in 
cavity QED [41]. 

Before we conclude, it should be mentioned that the idea of realizing multiqubit controlled phase gates with 
superconducting fiux qubits or charge qubits has been proposed previously [42,43]. Our present work, however, deals 
with the realization of a multiqubit controlled-C/ gate (a multiqubit controlled-" arbitrary transformation^^). Therefore, 
it is much more general than the previous works [42,43]. To the best of our knowledge, no one has yet demonstrated 
how to realize an n-qubit controlled- f7 gate within cavity QED. We believe that this work is of great importance since 
it provides a simple protocol to realize a multiqubit controlled-U' gate with SQUID qubits or atomic qubits within 
cavity QED. 
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Figure and Table Captions 



FIG. 1. Schematic circuit of an n-qubit controlled-?/ gate. A unitary transformation U is performed on the target 
qubit (qubit n) when the n — 1 controls on the filled circles (qubits 1, 2, and n — 1) are all in the state 
FIG. 2. Level diagram of a SQUID with the A-type three lowest levels |0) , |1) and |2). 

FIG. 3. (a) Resonant interaction of the cavity mode with the |0) ^ |2) transition of a SQUID, (b) Resonant 
interaction of the cavity mode with the |1) ^ |2) transition of the same SQUID. In (a) and (b), the coupling constant 
g is the same. Figure (b), where the ground level is |1) and the first excited level is |0) , is obtained by flipping Figure 
(a). The potential and the level structure shown in (a) and (b) can be obtained with external flux bias of (0.5 — 6) $o 
and (0.5 + 5) $0, respectively. 

FIG. 4. (a) Off-resonant interaction between the cavity mode and the |0) ^ |2) transition of a SQUID, (b) Off- 
resonant interaction between the cavity mode and the |1) ^ |2) transition of the same SQUID. The coupling constant 
g and the detuning A in (b) are the same as those in (a), which can be achieved with external flux bias of (0.5 — 6) $0 
and (0.5 + 5) $0 for (a) and (b), respectively. In (a) A = 1020 — uJc, while in (b) A = UJ21 — ijJc- 

FIG. 5. Change of the level structure (reduced) of SQUIDs (1, 2) during a two-qubit controUed-C/ gate performance. 
In (a), (b), (c), (d), (e), (f), (g), and (h), figures from left to right represent the level structures for SQUIDs 1 and 
2, respectively; the non-identical level spacings of the SQUIDs could be caused by nonuniform device parameters. In 
(a) and (c), the level spacings for the two SQUIDs are set to be much different, such that the irrelevant SQUID is 
decoupled from the applied pulse. The transition between any two levels linked by a dashed line is decoupled from 
the cavity mode, g and g are the off-resonant coupling constants between the cavity mode and the corresponding 
two-lc!vd transition of SQUID 2. gi and 52 arc the SQUID-cavity resonant coupling constants for SQUID 1 and SQUID 
2, respectively. The detuning A = u)2o — for (g) while UJ21 — oJc for (h). In addition, A = Wc — ^20- 

FIG. 6. Sketch of the setup for five SQUIDs (1, 2, 3, 4, 5) and a standing-wave quasi-one dimensional CPW cavity 
(Not drawn to scale). Each SQUID is placed in the plane of the resonator between the two lateral ground planes (i.e., 
the x-y plane) and at an antinode of the Be field. The two curved lines represent the standing-wave Be field, which 
is in the 2-direction. 

FIG. 7. Sketch of the setup for the realization of a controlled-t/ gate with three-level atoms and a cavity, (a) 
Resonant interaction of the cavity mode with the |0) ^ |2) transition of atom 1, atom I {I = 2, 3, n— 1), or atom n. 
(c) Resonant interaction of the cavity mode with the |1) <-^- |2) transition of atom n. (b), (d) Off-resonant interaction 
of the cavity mode with the |0) ^ |2) transition of atom n for different detuning setting, (e) Off-resonant interaction 
of the cavity mode with the |1) ^ |2) transition of atom n. gi, gi, and gn are the resonant coupling constants between 
the cavity mode and the |0) ^ |2) transition of atom 1, atom I {I = 2,3, ...,n — 1), and atom n, respectively, g'^ is 
the resonant coupling constant between the cavity mode and the |1) ^ |2) transition of atom n. and g!^ are 

the off-resonant coupling constantsjaetween the cavity mode and the corresponding two-level transition of atom n. 
A = W20 — 1^0 = ^21 — We, and A = Wc — W2o- 

TABLE 1. Parameters for a SQUID-cavity. f3L is the SQUID's potential shape parameter, R is the SQUID's 
effective damping resistance, and S is the surface bounded by the loop of the SQUID with width a and length b. ^2^ 
(7f^) is the energy relaxation time of the level |2) (|1)). 1^20 (i^2i) is the |0) ^ |2) (|1) ^ |2)) transition frequency. 

= {i\ $ |j) /$o is the magnetic dipole coupling matrix element between levels \i) and \ j) [i = 1, 2; j = 0, 1). I is 
the length of the quasi-one dimensional CPW cavity, A is the wavelength of the cavity mode with frequency fe, d is 
the gap between the center conductor and the adjacent ground plane, w is the width of the center conductor, t is the 
width of each ground plane, Lq is the inductance per unit length of the waveguide, and is the effective relative 
dielectric constant. 
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